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ABSTRACT 


Theoretical expressions are derived for the steady-state frequency response of 
a supported thermocouple wire. In particular, the effects of axial heat conduction 
are demonstrated for both a supported one material wire and a two material wire 
with unequal material properties across the junction. For the case of a one 
material supported wire, an exact solution is derived which compares favorably 
with an approximate expression that only matches temperatures at the support 
junction. Moreover, for the case of a two material supported wire, an analytical 
expression is derived that closely correlates numerical results. 

Experimental data were taken with a type K supported thermocouple. The 
test thermocouple was constructed with dimensions to demonstrate the effects of 
axial heat conduction assuming constant physical properties across the junction. 



L INTRODUCTION 


The evaluation of jet engine performance and fundamental studies of 
combustion phenomena depend on the measurement of turbulent fluctuating 
temperatures of the gas within the engine (Dils and Follansbee, 1976). Historically, 
these temperatures have been measured with thermocouples. The design of a 
thermocouple represents a compromise between accuracy, ruggedness and rapidity 

of response. 

For example, the measurement of fluctuating temperatures in the high speed 
exhaust of a gas turbine engine combustor is required to characterize the local gas 
density gradients or convective heat transfer (Fralick, 1985). Although 
thermocouples are suitable for the measurement of high frequency temperature 
fluctuations (< 1 KHz) in a flowing gas or liquid, the measured signal must be 
compensated since the frequency of the time dependent fluid temperature is 
normally much higher than the corner frequency of the thermocouple probe 
(Scadron and Warshawsky, 1952). Moreover, the amplitude and phase angle of the 
thermocouple response may be attenuated by axial heat conduction for rugged 
thermocouples of finite length (Elmore, et al.; 1983, 1986). 

In the present study, the theoretical steady-state frequency response of a 
supported thermocouple wire has been calculated to include the effects of axial 
heat conduction. These solutions, which represent an extension of earlier work 
(Forney and Fralick, 1991) are derived for both a supported thermocouple wire 
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with equal physical properties across the junction (e.g., roughly the same 
thermoconductivity, etc.) and a supported wire with unequal properties across the 
junction. Solutions are presented in the form of the amplitude ratio and phase 
angle for both cases. 

II. THEORY FOR ONE MATERIAL THERMOCOUPLE 

The steady-state frequency response of a thermocouple wire will be developed 
with the following assumptions: (a) the amplitude of the fluctuating fluid 

temperature is small relative to the mean absolute temperature (b) the 
thermocouple dimensions are small relative to the size of the turbulent eddies or 
enclosure dimensions (c) radial temperature gradients in a wire cross section can 
be neglected and (d) radiative heat transfer can be neglected relative to conduction 
and convection. 


In this section the geometry of Fig. 1 is considered where the material 
properties of thermal conductivity k, specific heat c and wire density pare assumed 
to be equal on both sides of the thermocouple junction. If the probe is immersed 
in a flowing fluid, the expression for the local conservation of energy in the 
thermocouple wire becomes (Scadron and Warshawsky, 1952) 


dT w d% Ah 
dx pcD 


at 


(T.-T.) 


( 1 ) 
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where a = k/ pc is the thermal diffusivity of the wire, T g is the ambient fluid 

temperature, h is the convective heat transfer coefficient, D is the wire diameter 
and T w is the local wire temperature measured along the axis at a distance x from 

the centerline (Fig. 1.) 

The wire and fluid temperatures are measured relative to the mean fluid 
temperature To. The ambient fluid temperature is taken to be a mean temperature 

together with a sinusoidal varying deviation from the mean, 

T t (t) = T 0 + T f e io ' (2) 

where co is the angular frequency of the ambient temperature. Since Eq. (1) is 
linear, we now seek a solution for the local wire temperature of the form 
(Hildebrand, 1976) 


T W =T 0 + f a {x)e i<a '. (3) 

Referencing all temperatures with respect to the mean gas temperature T 0 and 
normalizing with respect to the amplitude of the fluctuating ambient fluid 
temperature Tf, one defines a local normalized steady-state frequency response 

f(x) for the thermocouple wire of the form 


L-T 0 

T f 


= f( X y 


(4) 
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Substituting Eqs. (2) and (3) into Eq. (1), one obtains an ordinary differential 
equation of the form 


icoT 


= a 


d 2 T 

dx 2 


*<*>„(! -T) 


(5) 


where T = f(x) is the frequency response and the vector notation will be dropped 
for simplicity. Thus, for the geometry of fig. 1, one seeks a solution to the non- 
homogeneous linear second order differential equation for the dependent variable 
T of the form (Forney and Fralick, 1990) 


yT"~ G(co)T = -1. 


( 6 ) 


The general solution to Eq.(6) can be written in the form (Hildebrand, 1976) 


T(x) = 4sinh ^x + Bcosh ax + 

G(a) (7) 

where the parameters in Eqs. (5) and (6) are defined as 


co. = 


Ah 

pcD ’ 



G(w) = 1 + j' 


\°>nj 


(8) 


while in Eq. (7) the constants A and B are complex, 1 /G(co) represents the 
particular solution and the parameter 
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1. Approximate Solution 


Assuming that the material properties are constant across the junction and 
that the wire diameters are Di and D 2 in regions 1 and 2, respectively, Eq. (6) is 
subject to the boundary conditions 

TO = TO = T a (9a) 


T 2 (l + L ) = 0 . (9b) 

In this case, we seek a simple approximate solution that neglects the heat transfer 
at the interface between regions 1 and 2 at x = ± l where the parameters in Eqs. (5) 
and (6) are defined in terms of the wire diameters in each region. A similar 
approach will be used in a later section to obtain an approximate solution for the 
case in which the material properties of the two elements of the thermocouple are 
distinctly different. Hence, in region 1 





y = —, G^co) = 1+ i(— ) 
co, co, 


( 10 ) 


where co x is the natural frequency of the wire in region 1 of fig. 1. 


The solution to Eq. (6) for the one material wire on both sides of the junction 
in region 1 can be written in the form 

7](;c) = i4j sinh q x x + fi, cosh q x x + — - — . 

Ci (g>) (11) 


Substituting x = ±l in Eq. (11), the boundary conditions in Eq. (9) yield values for 
the constants 


A =0, 


B x = 


1 


cosh q x l 


T a 


n 

C,J‘ 


(12) 


Thus, one obtains a steady-state temperature distribution for the wire in region 1 
of Fig. 1 in the form 


t; W = 

Cj v cosh q x l ) ycoshqj ) (i 3 ) 

We now seek a solution in region 2 that satisfies the boundary conditions of 
Eq. (9). Since the temperature is symmetric about x = 0, it is convenient to define a 
continuous steady-state temperature distribution for the large wire of diameter D 2 
over the entire region ~(l+L) <x< ( l+L ) 0 r 


W = yr 
c 2 


1 cosh q 2 x 
cosh< 7 2 (/ + L) 


( 14 ) 
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Since T a = T 2 (l), one obtains the boundary value from Eq. (14) 


T = — 

° g 2 


_ cosh< 7 2 / 
cosh< 7 2 (/ + L) 




( 15 ) 


Substituting the value for T, at x = l of Eq. (15) into Eq. (13), the approximate 
temperature distribution in region 1 becomes 



f cosh q x x y 

1 f 
+ 

K cosh qj j 

G 2\ 


cosh< 7 2 / 


cosh< 7 2 (/ + L) J cosh qj 


cosh<7iX 


(16) 


Thus, the approximate frequency response at the thermocouple junction (x = 0) for 
the one material wire becomes 


VO) = -i- 

G l 


1 - 


1 ] 

1 

i 

r cosh q 2 l ^ 

1 

cosh qj j 

c 2 

< cosh q 2 (l + L)j 

cosh q x l 


(17) 


The steady-state frequency at the thermocouple junction x=0 is normally 
characterized graphically in the form 


J 


T(0) = \T(0)\e‘* 


(18) 
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where 1^(0)! is the amplitude ratio and <j> is the phase angle. In the latter case, the 
phase angle in degrees is 


<J> = 57.3 tan 1 



(19) 


where Im[T(0)] and Re[T(0)] are the imaginary and real parts of T(0), respectively. 

2. Exact Solution 

If the boundary conditions listed in Eq.(9) include equal rates of conductive 
heat transfer at the interface between the thermocouple and support wires at 
x = ± l, the exact solution is subject to 


71(0 = T 2 (I) 


(20a) 


1 dr 2 dx 


(20b) 


T 2 (l + L) = 0 . 

Since the solution to Eq. (6) in region 1 is of the form 


(20c) 


T x (x) = sinh q x x + B i coshq l x + — 

G i 


( 21 ) 
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where by symmetry ^I(0 _ ^I(“0/ one obtains Ai - 0. Thus, the form of the 
solution in region 1 becomes 


^(x) = Bjcosh^x + — . 


( 22 ) 


For region 2, where the spatial coordinate is in the range l£x£(l + L) f it is 
convenient to write 


T 2 (x) = A 2 s\nhq 2 (I + L-x) + B 2 cosh q 2 (l + L-x) + 



(23) 


From the boundary condition T 2 (l + L) — 0, one obtains B 2 — l/G 2 or the form of 
the solution in region 2 becomes 

T 2 (x) = A 2 s\nhq 2 (l + L-x) + ^-[l - cosh <? 2 (/ + L-x)] . ^ 


Substituting Eqs. (22) and (23) into boundary condition (20a), one obtains a 
linear equation for the constants B 1 , \ or 


cosher,/ 


i4jSinh<7 2 L = 


-^-[1-cosh^L] - -U 

u 2 t/, 


(25) 


Similarly, substituting Eqs. (22) and (23) into boundary condition (20b), one obtains 
a second linear equation for B v \ or 
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2?,2sinh<7,/ + cosh q 2 L - — sinh q 2 L 

G 2 


( 26 ) 


where the complex constant Q is defined as 


D 2 2 q 2 


(27) 


Solving Eqs. (25) and (26) for Bi and A 2 , one obtains the determinate system 


4 = 


J-[l-cosh % L]-i- 
C/ 2 G x 

— sinh<7 2 L 

G 2 


-sinh<7 2 L| 
cosh q 2 L 


DET 


(28) 


and 


*z = 


cosh <7,1 
\Q sinh qj 


^-[l-cosh^LJ-i-l 


-5-sinh<7 2 L 

g 2 


G, 


DET 


(29) 


where the determinate in the denominator is equal to 


DET = coshq,/ cosh<7 2 L + Qsinh^/sinh^L. 


(30) 


Thus, solving for the constants Bi and A 2 from Eqs. (28) and (29) and substituting 
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Bi into Eq. (22), one obtains an exact expression for the steady-state temperature 
distribution in the form 


= 


— (cosh<7 2 L-l) - -^-cosh<? 2 L 
G 2 G l 


cosh <7^ 


cosh <7^ cosh q 2 L + Qsinh qj sinh q 2 L 


1 

— + — . 
G x 


(31) 


Thus, the steady-state frequency response at the thermocouple junction x=0 
becomes 


7(0) = 


-^-(cosh^ 2 L- 1) - -^-cosh^ 2 L 


cosh <7^ cosh q 2 L + <2 sinh q x l sinh q 2 L 



(32) 


IE. THEORY FOR TWO MATERIAL THERMOCOUPLE 

Certain types of thermocouples have distinctly different material properties 
across the junction. For example, a copper-constantan thermocouple has a 
thermal conductivity on the copper side that is more than an order of magnitude 
larger than constantan. In this case, the expressions developed in the previous 
section for the frequency response are in error since unequal material properties 
would provide an asymmetric temperature profile. 

Referring to fig. 2, the thermocouple schematic now has four distinct regions 
that are distinguished by either different wire diameters or physical properties. For 
example, on the left side of the schematic of fig. 1 the thermal conductivity, density 
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and specific heat have the values ki. Pi and ci , respectively, while on the right side 
of the schematic the material properties are V. 2 , p 2 and C 2 - 

1. Temperature Distribution for Small Wire 

Since the differential equation describing the steady-state frequency response 
Eq. (6) applies to all regions of the schematic of fig. 1, the steady-state frequency 
response in regions 1 and 2 are given, respectively, by the expressions 

T x (x) = /ijSinh^x + A cosh <7,* + — 

G i (33) 

T 2 (x) = sinh q 2 x + B 2 cosh q 2 x + — . 

G i (34) 

Here, the four constants designated by A and B in Eqs. (33) and (34) are determined 
by the four boundary conditions 

OT = T 2 ( 0) (35a) 


up 2 = I’D 2 ^ 2 ^^) 

** 1 dx ** 1 dx 


(35b) 


r iH) = r a (35c) 

r 2(0= T b • (35d) 
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Solving for the four values of the constants designated by A and B in Eqs. (33) 
and (34), one obtains an expression for the steady-state temperature distribution in 

region 1 of the form 



|sinh^,(x+/) 


, 0 * OJ 


osh<? 2 / + | T b 


(’ 



+ 



A 



\) 


sinh^xcosh^/ 



(36a) 


and in region 2 


T(x) = jr + [~Q t sinh^(/- 


*) 


IK G 2 G \) 


1 ^ 

z°^ql- \T q ~ — 
1/ J 


( 


1 




T 

6 G 2 y 


^cosh<7 2 x sinh^/ + Q { cosh ql sinh <7 2 *j| 




(36b) 


Here, it should be noted that the two functions above can be obtained from each 
other by interchanging the subscripts 1 and 2, the values a and b and by changing 
the sign of x. Also, in Eqs. (36a) and (36b) the parameters 


Aj = Q, cosh <7,/ sinh q 2 l + cosh q 2 l sinh q x l 


2. Temperature Distribution for Large Wire 

In this section a solution is sought for the large wire in regions 3 and 4 of the 
schematic of fig. 2. To simplify the analysis, a continuous solution is derived for 
the steady-state frequency response over the entire region -(/ + £)£*<;(/+£,). This 
solution must satisfy the boundary conditions 

(37a) 


1 3 dx 2 3 dx 


(37b) 


r(-/-L) = o 


(37c) 


r 4 ( / + L )-°- ( 3 7d) 

In principle, it is now possible to derive a continuous solution for T$(x) over 
the range -{1 + L)<x < 0 and for T4M over the range 0 < x < (/ + L). These solutions 
are obtained from Eqs. (36a) and (36b) by replacing the subscripts 1 — » 3 and 2 — > 4 
and by redefining the parameters /->/ + £ and setting the boundary conditions 
T * =T b = 0- Moreover, the boundary condition 7’ a =T 3 (-/) i s obtained from the 

resulting expression for ^(x) by substituting x = —l while the boundary condition 
Tb = T4O). Thus, one obtains values for T a , T b that appear in Eqs. (36a) and (36b) in 
the form 
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( 38 ) 


t = i-+{sinh q 3 L 


LV G 4 C 3 J 


cosh^ 4 (/ + L)--~ 

^4 


- —\Q, cosh <? 3 / sinh <7 4 (/ + L) + cosh q A (/ + L) sinh * 3 ^])^ ^ 


and 


T b = ”7* + [Q t sinh q 4 L 

^4 


---1 

, C 3 C 4 J 


cosh<7 3 (/ + L)--^- 


_ -L [cosh <7 4 / sinh <7 3 (/ + L) + Q, cosh<7 3 (/ + L)sinh<7 4 /]}[ -j- 
G 4 l v 2 


(39) 


Here, the parameters are defined as 


and 


A, = Cj cosh q i (/ + L) sinh q A (/ + L) + cosh q 4 (l + L) sinh <? 3 (/ + L) 


a - p- 


3. Frequency Response 

The steady-state frequency response for the two material supported 
thermocouple is obtained from Eq. (36a) or (36b) by setting x — 0. 
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no) 


= sinh< 7 2 / r - + ^-(cosh< 7 j/-l)J + sinh^jV*-!- ~ (cosh^ 2 / - l) 


(40) 


where T a and Tb are given by Eqs. (38) and (39). 


It should be noted that the steady-state frequency response Eq. (40) provides a 
wire temperature that is continuous everywhere and conserves the heat flux at the 
junction x = 0. This represents an approximate solution since the heat flux at the 
interface between the large and small wires x = ±l has been neglected. 

IV. RESULTS 


The amplitude ratio and phase angle of the thermocouple frequency response 
were plotted graphically for the case of a one material wire as shown in the 
schematic of fig. 1. In this case, average properties of a type B or 
Pt / 6% Rh-Pt / 30% Rh were used since the material properties were nearly equal 
across the thermocouple junction. The wire dimensions, properties and gas 
conditions are listed in table l(Touloukian et al., 1970). 

The amplitude ratio and phase angle were also plotted for a two material 
thermocouple wire as shown in the schematic of fig. 2. In this case, a type B 
thermocouple described in table 1 was used in addition to a type T or copper - 
constantan described in table 2 (Touloukian et al., 1970). 

The form of the convective heat transfer coefficient h that appears in the 
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computation of the natural frequency defined in Eq. (8) was determined from 
the expression (Scadron and Warshawsky 1952) 

Nu = .485Re^Pr^ (41) 

where Nu(=hD/k f ) is the Nusselt number, kf is the thermal conductivity of the 

ambient fluid, Pr(= *>//«) j s t he Prandl number, and Re (=vD/v f ) j s the 

Reynolds number of the thermocouple wire. Here, v and v f are the fluid velocity 
and kinematic viscosity, respectively. It should be noted that the convective heat 
transfer coefficient h a D m and the natural frequency of a thermocouple wire for 
given material properties &>* a D m • 

1. One Material Thermocouple 

The amplitude ratio |7X0)| at the wire junction for the steady-state frequency 
response derived from Eqs. (17) and (32) is shown in fig. 3. This assumes a type B 
thermocouple wire with the dimensions listed in table 1. In this case, the average 
material properties listed in table 1 were used since the one material theory 
assumes that the properties of the thermocouple wire are uniform across the 
junction. It is evident in fig. 3 that the amplitude ratio derived from the 
approximate expression Eq. (17) is nearly identical to the exact derivation Eq. (32). 
Thus, it appears that the conservation of heat flux at the interface between the 
small and large wires of the schematic of fig. 1 is of secondary importance. 

Also included in fig. (3) are the numerical computations of Stocks (1986). 
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These solutions represent explicit finite difference solutions to the one 
dimensional unsteady heat transfer equation as shown in Eq. (1). The small 
deviation of the numerical results from the exact solution at low frequency in fig. 
(3) is apparently due to the unsteady character of the numerical results. Similar 
computations of the phase angle $ for the type B thermocouple are represented in 
fig. 4. As indicated, the phase angle varies over the range 0£<J>£-7r/2 and 
approaches the lower limit of -Jt/2 for large frequencies <*>/ o* » 1 . 

The spacial variation of the amplitude ratio |7\*)| derived from Eqs. (17) and 
(32) is graphed in fig. 5. These computations were made at an angular frequency of 
= 0.1 for the type B thermocouple. As evident in fig. 5, the difference 
between the exact and approximate expression is somewhat exaggerated at a very 
low frequency. Nevertheless, the error represented by the approximate solution is 
less than 7% over the length of the thermocouple. As stated earlier, matching the 
heat flux at the interface between the small and large wire at jc = ±/ in the 
schematic of fig. 1 appears to be of secondary importance in relation to providing a 
continuous temperature profile along the wire. 


The amplitude ratio |^(0)|, phase angle <j> and spacial variation Y are also 
plotted in figs. (6), (7) and (8), respectively, from the steady-state frequency response 
represented by Eqs. (17) and (32). These results represent a type B thermocouple 
with a smaller diameter ratio D 2 /D l = 2 (see table 1 for dimensions). As indicated 
in figs. (6), (7) and (8), the approximate and exact solutions represent comparable 
results in all three cases. Thus, the diameter ratio of the large and small wire 
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appears to have little effect in the favorable comparison between the approximate 
and exact solutions representing the one material steady-state frequency response. 


2. Two Material Thermocouple 

The amplitude ratio | r H at the wire junction for the steady-state frequency 
response derived from Eqs. (17) and (40) is shown in fig. 9. This assumes a type B 
thermocouple wire with the dimensions listed in table 2. In this case, the average 
material properties listed in Table 2 were used for the amplitude of the frequency 
response derived from the one material solution of Eq. (17). Also plotted in fig. 9 
is the amplitude ratio derived from the two material solution of Eq. (40). In the 
latter case, the individual material properties listed in Table 2 were used. 

As expected, the amplitude ratio for the steady-state frequency response of a 
type B thermocouple is nearly identical with either the one material or two 
material approximate solutions. This is a consequence of roughly equal material 
properties across the junction for type B thermocouples. This plot also validates 
the two material approximate solution Eq. (40) and the values of the boundary 
conditions for T a and T b substituted from Eqs (38) and (39). The same conclusion 
can be drawn with respect to the phase angle <J>, shown in fig. 10. 

The spadal variation of the amplitude ratio |7T*)| derived from Eqs. (17) and 
(40) is plotted in fig. 11. These computations were made at an angular frequency of 
(of(o n = 0.1 for the type B thermocouple with the dimensions listed in table 2. It is 
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interesting to note the small asymmetry in the amplitude ratio on the left and 
right side of the junction. This asymmetry is the result of the small differences in 
the physical properties across the junction as listed in table 2. 

The steady state amplitude ratio F(0)| for a type T thermocouple is plotted in 
fig. 12 using equations (17) and (40). The dimensions and material properties are 
listed in table 2. The average material properties listed in Table 2 were used to 
compute the amplitude ratio of the frequency response with the one material 
solution, Eq. (17). Also plotted in fig. 12 is the amplitude ratio derived from the 
two material solution, Eq. (40). In the latter case, the individual material 
properties also listed in table 2 were used. 

As indicated in fig. 12, the amplitude ratio for the approximate one material 
steady-state frequency response of a type T thermocouple is distinctly different 
from the two material approximate solution. This is a consequence of unequal 
material properties across the junction for the type T thermocouple. Also shown 
in fig. 12 is a numerical solution of the second order ordinary differential equation 
for the temperature, Eq. (6). The numerical finite difference solution of the 
boundary value problem of Eq. (6) matches both the temperature and heat flux at x 
= 0 and * = ± / in the schematic of fig. 2. It is clear from fig. 12 that the approximate 
two material analytical solution, equation (40), is in good agreement with the 
numerical results despite some differences at low frequencies (oJa) x < .04. Similar 
results are indicated for the phase angle <i>in fig. 13. 

The spadal variation of the amplitude ratio iT'fx)! derived from Eqs. (17) and 
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(40) is graphed in fig. 14. These computations were made at an angular frequency 
of to/cOy =0.1 for the type T thermocouple. As evident in fig. 14, the two material 
solution derived from Eq. (40) accurately represents the features of the asymmetry 
associated with a type T thermocouple. In particular, the relatively large resistance 
to axial heat conduction in the constantan wire on the right of the junction is 

reflected in the larger values of the amplitude ratio PX-*)|. Also shown in fig. 14 is 
the approximate one material solution represented by Eq. (17). In the latter case, 
the average values for the material properties of a type T thermocouple were used 
as listed in Table 2. Therefore, one can conclude that for thermocouples whose 
material properties on either side of the junction are markedly different, the two 
material solution developed in this paper is a substantial improvement in 
accuracy both in frequency response and in temperature distribution along the 
wire. 

V. Experimental Procedure 

In the experiment described below, thermocouple sensors were exposed to a 
constant velocity air stream of varying temperature. In particular, the dynamic 
response of the thermocouple sensors to a square wave temperature profile was 
measured for a range of frequencies. 

1. Rotating Wheel Experiment 

A rotating wheel configuration was used to deliver the test airstream to the 
proposed sensors. A similar experimental apparatus is described in detail by 
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Elmore et al. (1986). A schematic of the rotating wheel apparatus used in the 
present experiment is shown in figure 15. As the wheel rotates, the holes pass the 
two air supply tubes which allow slugs of hot and cold air to alternately enter a 
transition tube assembly mounted directly above the rotating wheel. In the 
transition tube the slugs of hot and cold air coalesce into a single air stream 
providing roughly a square wave temperature profile covering a range of 
frequencies from roughly 1 to 30 Hz. 

The analog temperature signal was digitized with a Data Translation DT-2801 
board mounted in an expansion slot of an IBM AT compatible computer as shown 
in figure 15. The ASYST software loaded on the hard disc of the personal 
computer provided a flexible system for data storage, manipulation and display. 

The true temperature profile of the airstream was measured with a constant 
current anemometer (TSI 1054- A) and sensor (1226 PI 2.5). 

2. Thermocouple Construction 

Thermocouple wire of the desired length and type is threaded into the four 
hole ceramic with the thermocouple end last. Three or four kinks are made in 
each wire near the thermocouple end so the wire must be firmly pulled into the 
tube leaving enough wire sticking out to make the thermocouple. Drops of epoxy 
are picked up with a piece of .010" diameter wire, added to the ceramic tube at A 
(see fig. 16) and pushed down around both wires. The kinks and epoxy firmly 
fasten the wires in the ceramic so they do not twist when the free ends are 
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manipulated for electrical connections. 


The junction at B in figure 16 is made by cutting the large wire about half way 
through with a razor blade, laying the small diameter wire in the cut, and welding 
the cut closed with a stored energy spot welder. After the weld, the excess wire 
ends that protrude through the weld area are bent and broken off to clean up the 
junction. 

For the chromel alumel couples, the junction at C of figure 16 was made using 
a stored energy spot welder. The wires are crossed and welded and then trimmed 
with a razor blade and bent with tweezers under a microscope until collinear. 

A stored energy spot welder will not work for the copper constantan couples. 
For these couples the junctions were made by silver soldering. The silver solder 
wire is coated with flux and the end heated with a torch until a drop forms. A 
twisted pair of thermocouple wires is pushed briefly into the flux coated solder 
drop. The solder will wet the pair up to where the twist stops. Again using a razor 
blade and tweezers under a microscope, the thermocouple is bent and trimmed 
until collinear. 

3. ASYST Software 

ASYST software was developed to acquire temperature data simultaneously 
from the thermocouple and constant current anemometer. Data were digitized for 
two channels at a sampling rate of 512 Hz per channel for a total sample time of 
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two seconds. The ASYST software code is listed in figure 17. 

The ASYST code first plotted the temperature profiles versus time from the 
thermocouple and anemometer and an example of the plot is shown in figure 18. 
The ASYST software next took the Fast Fourier Transform of the temperature data 
in each channel and recorded the amplitude ratio and phase angle between both 

channels at the first harmonic for the square wave. These data are discussed in the 
next section. 

VI. Results 

Initial tests with a signal generator and the output from an RC circuit indicated 
that the ASYST code and data acquisition hardware were operating properly. 
Several test sensors were constructed and tested with the dimensions listed in 
table 3. In each case measurements of the amplitude and phase angle were 
compared with the appropriate theory. A discussion of the results is given below. 

1. First Order Response 

The amplitude ratio and phase angle were measured with the type K 
thermocouple listed first in table 3. The lengths of the thermocouple and support 
wires in this case were chosen to eliminate the effects of axial heat conduction. 
The experimental data representing the amplitude ratio are plotted in figure 19 
along with the theory representing a first order frequency response. The 

experimental data were correlated with a natural frequency <y, =5.5 sec" 1 and a 
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thermal diffusivity of cc =.059 cm 2 / sec. 


As indicated in figure 19 the data has been correlated with the theory to within 
10%. Although the experimental data is reproducible, there is a noticeable drift of 
the data relative to the theory with increasing wheel frequency to- The drift is 
probably due to a small increase in the natural frequency ©j of the thermocouple 
wire with increasing wheel frequency a. 

The natural frequency of Equation (10) depends on accurate predictions of 
the convective heat transfer coefficient h from Equation (41). However, there may 
be a small change in either the mean air velocity or unsteadiness in the mean 
velocity which would lead to small changes in the heat transfer coefficient h and 
the natural frequency (0 ,. 

Also plotted in figure 20 are experimental data representing the phase angle <t> 
for the first order frequency response of the type K thermocouple. In this case the 
data is correlated to within 7% of the theory. The experimental data for the phase 
angle also indicate a small drift in the natural frequency as was apparent in 
figure 19. 

2. Frequency Response (one material) 

The amplitude ratio and phase angle were measured with the second type K 
thermocouple listed in table 3. The lengths of the thermocouple and support 
wires in this case were chosen (i.e., shortened) to demonstrate the effects of axial 
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heat conduction. The experimental data representing the amplitude ratio are 
plotted in figure 21 along with the theoretical prediction given by Eq. 32. Both the 
data and theory were correlated by assuming that the material properties are 
uniform across the thermocouple junction. In this case the data were correlated 
with a natural frequency eo } = 5.0 sec -1 and a thermal diffusivity of 
a - .059 cm 2 / sec. 

As indicated in figure 21 the data has been correlated with the theory to within 
10% over most of the range of wheel frequencies. The data is reproducible but 
again a noticeable drift of the experimental data exists relative to the theory as was 
discussed in the previous section. 


Also plotted in figure 22 are experimental data representing the phase angle <j> 
for the frequency response of the shortened type K thermocouple. The data is 
correlated to within 7% of the theory. The experimental data for the phase angle 
also indicate a small drift in the natural frequency ,. 
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IV. 


NOMENCLATURE 


A 

B 

c 

D 

G 

h 

i 

k 

kf 

/ 

L 

Nu 

Pr 

Q 

Qt 

Qs 

q 

Re 

t 


constant of integration 
constant of integration 

material specific heat [J-gm x - °^~ 1 
thermocouple wire diameter (cm) 

1 +i(®/®„) 

heat transfer coefficient (/-cm 2 ~ 5 


unit imaginary number ( = 

material thermoconductivity (/-cm 1 -s 1 - K ) 

gas thermoconductivity (/-cm 1 - s 1 - °K ’) 
length of small thermocouple wire (cm) 
length of large thermocouple wire (cm) 

Nusselt number (= hD/kf) 

Prandl number ( = v // a ) 

£& 

oli h 


fyt 2 
K Si 

(G/y ) 5 

Reynolds number (= vD/u f ) 
time (s) 
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T 


steady-state frequency response 

T f 

= 

amplitude of periodic gas temperature (°^) 

T g 

— 

gas temperature (°^) 

To 

= 

mean gas temperature (°^) 

—* 

T 

A a> 

= 

complex amplitude of periodic wire temperature 1 

T 

1 w 

= 

local wire temperature (°^) 

V 

= 

gas velocity (cm -s' 1 ) 

X 

=: 

axial distance from center of wire (cm) 

Greek Symbols 


a 

= 

thermal diffusivity ( cm 2 -s~ l ) 

Y 

= 

a/co„ (cm 2 ) 

v r 

= 

kinematic viscosity of gas ( cm 2 ~ fl ) 

CO 

= 

angular frequency ( s ') 


= 

natural frequency of wire (= 4h/pcD)(s~') 


= 

phase angle 

P 

= 

.3 

material density (gm - cm ) 

\ 

= 

dimensionless function 

a 2 

= 

dimensionless function 
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Table 1 - Properties of One Material Wire (Type B) 


5i 

.025 

.0076* 

Average Properties 



3.8 


*<»! = 32.9 sec' 1 


Dimensions (cm) 


52 

.05 

.038 


52^5l I L 

2 .2 .35 

5 .1 .2 


F Type B 
r cm 2 > | 

.sec J 


Gas Properties 


T 0 = 900 °AT 


.22 


M= .26 
P = 1 atm 
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Table 2 - Properties of Two Material Wire 


Dimensions (cm) 

Dj Pz/Pl- i L 

.0076 .038 5 .1 .2 


Air Properties 
T 0 = 300 °K 
P= 1 atm 
V = 50 m / sec 


Properties of Type B 



P {cm 3 -o*] 

4—1 

|_ sec J 

Pt - 6% Rh 

2.73 

.238 

Pt - 30% Rh 

2.86 

.190 

Average 

2.8 

.214 


Properties of Type T 



Copper 

3.44 

1.16 

Constantan 

3.48 

.067 

Average 

3.46 

.614 


Wire Location 


Region 

Type B 

Tvpe T 

1 

Pt - 6% Rh 

Copper 

2 

Pt - 30% Rh 

Constantan 
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Table 3 - Dimensions for Test Thermocouples (cm) 


Test 

Thermocouple 

Type 

Eh. 

Efc 

Ug/Di 

i 

L 

First Order 

K 

.0076 

.038 

5 

0.25 

0.5 

Axial Conduction 
(one material) 

K 

.0076 

.038 

5 

0.1 

0.2 
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Fluid Flow 


1- 0 
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Type B (D 2 /D 1 = 5) 
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le rati 



Type B (DjD, = 2) 
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Type B (D 3 /D 1 = 5) 
— two material 
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Amplitude ratio vs. angular frequency for two material wire. Solid line 
is Eq.. (40). Crosses are approximate one material solution Eq. (17). 
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11. Amplitude ratio vs. distance for two material wire. Solid line derived 

from Eqs. (36a) and (36b). Crosses are approximate one material solution 
Eqs. (14) and (16). 



12. Amplitude ratio vs. angular frequency for two material wire. Solid line 
is Eq. (40). Crosses are approximate one material solution Eq. (17). 
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13. Phase angle vs.angular frequency for two material wire. Solid line is Eq. 
(40). Crosses are approximate on material solution Eq. (17). 
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HIGH PRESSURE 
AIR 



15. Rotating wheel experiment. 
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c 



16 . Schematic of test thermocouple 
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FORGET. ALL 

REAL DIME 1024 ] ARRAY TIME 

COMPLEX DIME 1024 ] ARRAY TRANS 

REAL DIME 6 ] ARRAY ANS 

REAL DIME 1024 ] ARRAY ZMAGO 

REAL DIME 1024 ] ARRAY ZARGO 

INTEGER SCALAR NUM 

REAL SCALAR MEN 

REAL SCALAR CD 

REAL DIME 1024 ] ARRAY FREQS 

REAL DIME 1024 ] ARRAY SIGNAL 

INTEGER dimE 1024 , 2 ] array DATA. BUFFER 


1. CD := 

.0000 DATA. BUFFER := 

LOAD. OVERLAY ACQUIS. SOV 
.0000 1.0000 A/D. TEMPLATE 

DATA. BUFFER TEMPLATE . BUFFER 
CYCLIC 

2048.0000 TEMPLATE . REPEAT 
CD CON VERS ION. DELAY 
DEMO. TEMPLATE A/D.INIT 
A/D. IN>ARRAY 
LOAD. OVERLAY WAVEOPS.SOV 
1024 REAL RAMP 1 - 2 * CD * 


DEMO . TEMPLATE \ DEFINE AN A/D TEMPLATE 

\ DECLARE ARRAY AS A TEMPLATE BUFFER 
\ SET TEMPLATE BUFFER TO CYCLIC MODE 
\ SET REPETITIONS FOR I/O INSTRUCTION 
\ CONVERSION RATE (MSECS/SAMPLE) 

\ INITIALIZE CURRENT A/D TEMPLATE 
\ A/D INPUT TO TEMPLATES BUFFER 

1000 / TIME := \ SET TIME AXIS 


DATA. BUFFER XSECT[ 1 ] 2048 - \ 

SIGNAL := SIGNAL 

-2048 / 10 * -55.274 * 353. + \ 

SIGNAL := SIGNAL 

0 . 5 SET . CUTOFF . FREQ \ 

SMOOTH SIGNAL := SIGNAL 
MEAN MEN := SIGNAL 

MEN - SIGNAL := \ 

TIME SUB[ 1 , 500 , 1 ] SIGNAL SUB[ 
XY . AUTO . PLOT S I GN AL \ 

FFT \ 

TRANS := TRANS ZMAG ZMAGO := \ 


CHANNEL 1 ON STACK (ANEMOMETER) 
CALIBRATE (DEGREES C) 

SMOOTH DATA (CYCLES/POINT) 


CENTER ON ORIGIN 
1 , 500 , 1 ] 
PLOT CHANNEL 1 
TAKE FFT 

MAGNITUDE OF FFT 


5 SET.#. POINTS 
1 SET.#. OPTIMA 

ZMAGO SUBC 1 , 120 , 1 ] LOCAL. MAXIMA \ FIND INDEX AND MAX OF MAGNITUDE 
SWAP NUM •= NUM ' \ INDEX OF MAX MAGNITUDE 

1 - CD / 1024. / 1000 * PI * \ FREQUENCY AT MAX MAGNITUDE 

ANS [ 1 ] := ANS [ 1 ] DROP ANS [ 2 ] := 

TRANS [ NUM ] DEG ZARG ANS [ 3 ] := \ ARGUMENT AT MAX MAGNITUDE 


DATA. BUFFER XSECT[ 2 ] 2048 - \ CHANNEL 2 ON STACK (THERMOCOUPLE) 

SIGNAL := SIGNAL 

.124 * SIGNAL := SIGNAL \ CALIBRATE (DEGREES F) 

0.5 SET. CUTOFF. FREQ \ SMOOTH DATA (CYCLES/POINT) 

SMOOTH SIGNAL := SIGNAL 
MEAN MEN := SIGNAL 


11:42:33.77 


ASYST Version 3.00 

Page 1 TEMP3 . FOR 03/13/91 
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MEN - SIGNAL := 

TIME SUBC 1 , 500 , 1 ] 

SIGNAL SUB[ 1 , 500 , 1 ] 
XY. DATA. PLOT SIGNAL 
FFT 

TRANS := TRANS ZMAG ZMAGO 
5 SET.#. POINTS 
1 SET.#. OPTIMA 


\ CENTER ON ORIGIN 


\ PLOT CHANNEL 2 
\ TAKE FFT 
\ MAGNITUDE OF FFT 


SUSC 1 , 120 , 1 ] LOCAL. MAXIMA \ FIND INDEX AND MAX OF MAGNITU 
SWAP NUM := NUM \ INDEX OF MAX MAGNITUDE 

J “ ? D / 1024. / 1000 * PI * \ FREQUENCY AT MAX MAGNITUDE 

ANS C 4 ] := ANS [ 4 ] DROP ANS [ 5 ] := 

TRANS [ NUM ] DEG ZARG ANS [ 6 ] := \ ARGUMENT AT MAX MAGNITUDE 


CR 

FREQ MAG PHI " CR 

ANS [ 1 ] ? CR 

ANS [ 5 ] ANS [ 2 ] / ? CR 

ANS C 3 ] -1 * ANS [ 6 ] + 180 PI / CD * 1000 / ANS [ 1 ] * - ’ CR 
FORGET. ALL 


ASYST Version 3.00 

Page 2 TEMP3. FOR 03/13/91 11:42:34.65 


17. ASYST software 
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.059 cm 
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20. Phase angle vs. angular frequency. ^ = 5.5 sec \ a- .059 cm 




54 


ORIGINAL PAGE IS 

OF POO* QUALITY 


21. Amplitude ratio vs. angular frequency. Solid line is Eq. (32) 

m ^ n ~ — ncn 2 y * x ' 
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Phase angle vs. angular frequency. Solid line is derived from Eq. (32). 
ft), = 5.0 sec -1 , a = .059 cm 2 / sec. 





